Live Soap: Order, Fluctuations and Instabilities in Active Smectics 



Tapan Chandra Adhyapak, 1, * Sriram Ramaswamy, 1, t and John Toner 2 ' * 

1 Centre for Condensed Matter Theory, Department of Physics, 
Indian Institute of Science, Bangalore 560 012 India 
2 Department of Physics and Institute of Theoretical Science, University of Oregon, Eugene, OR 97403, USA 

(Dated: March 2, 2013) 

We construct a hydrodynamic theory of noisy, apolar active smectics, in bulk suspension or on a 
substrate. Our predictions include: quasi-long-ranged smectic order in dimension d = 2, and long- 
ranged in d = 3, extending previously published results to all dynamical regimes; Kosterlitz-Thouless 
melting to an active nematic at high and low concentrations in d — 2; nonzero second-sound speed 
parallel to the layers; the suppression of giant number fluctuations by smectic elasticity; instability 
to spontaneous undulation and flow in bulk contractile smectics; a layer spacing instability, possibly 
oscillatory, for large enough extensile active stresses. 



Active particles [1] in a state of orientational order 
[2, 3] exhibit fluctuations [4-7] and flow properties [8-13] 
that differ strikingly from those in equilibrium systems 
with the same spatial symmetry. Translationally ordered 
active matter has received less attention [10, 14]. This 
paper studies active systems with spontaneously broken 
translation-invariance in one direction - active smectics 
- in a broad range of dynamical regimes. 

We consider apolar systems of particles carrying an 
axis of orientation disposed on average along the normal 
to the smectic layers (i.e., Smectics A [15]), with active 
stresses [11] pulling in or pushing out, i.e., contractile or 
extensile, along that axis. A variety of such models with 
this symmetry are possible, depending on conservation 
laws and geometry; we study five of these. The simplest 
is systems with no conserved quantities. Such a model 
describes the dynamics of Rayleigh-Benard roll patterns 
[16, 17], and spontaneously layered flocks of self-propelled 
apolar entities, reproducing or dying while in motion [18], 
on a substrate which serves as a momentum sink. The 
second is layered flocks moving on a substrate with num- 
ber conservation, but without momentum conservation. 
The third is bulk layered systems in a background fluid 
with momentum conservation treated in the "Stokesian" , 
i.e., viscosity-dominated, limit appropriate for colloidal 
or microbial active systems, and the fourth is such sys- 
tems confined between no-slip walls, where the surfaces 
are a momentum sink and the hydrodynamic interaction 
is screened at long wavelengths. We conclude by analyz- 
ing bulk systems in a background fluid at wavelengths 
beyond the Stokesian regime, where inertia dominates 
over viscosity. 

Our results: (i) In all the five cases we study, smec- 
tic order, when dynamically stable, is long-ranged in the 
presence of noise in dimension d — 3 and quasi-long- 
ranged in d = 2. This reinforces and extends the findings 
of [14]. (ii) The active smectic undergoes a transition to 
an active nematic as the concentration of active parti- 
cles is varied, in all five cases. In two dimensions, "reen- 
trance" [19] necessarily occurs: the active nematic occurs 
at both large and small concentration, with the active 



smectic at intermediate concentrations. Both transitions, 
in the "no- conservation" case, are of Kosterlitz-Thouless 
type [20]; the nature of the transitions for the other four 
models is unknown, (iii) Bulk smectic liquid crystals in 
the Stokesian limit are hydrodynamically stable to the 
presence of extensile active stresses with magnitude be- 
low a threshold value [21] (iv) Active smectics, unlike 
their oricntationally ordered counterparts [2-4, 6], have 
finite concentration fluctuations, whose magnitude, how- 
ever, diverges, as activity approaches the threshold, (v) 
Bulk active smectics with contractile active stresses are 
generically unstable without threshold to spontaneous 
undulations [22] and flow, (vi) Confined active smec- 
tic suspensions are in general stable at long wavelengths 
for small enough active stresses of either sign. Beyond 
a threshold value of activity the confined system too 
undergoes an instability, which is likely to be oscilla- 
tory for the extensile case, (vii) For stable bulk active 
smectics, the speed of the smectic second sound mode 
is nonzero for propagation parallel to the smectic layers. 
Our findings apply to active smectics in a wide range 
of settings including vibrated granular layers [23] and 
the Rayleigh-Benard problem [16]. Agitated 2DEGs [24], 
where Coulomb and magnetic-field effects enter, will be 
discussed elsewhere [25]. 

We begin with the simplest case: active elements whose 
number is not conserved, spontaneously condensed into a 
uni-directional periodic structure, i.e., a smectic A, with 
mean layer normal n along z. The only hydrodynamic 
field in this case is the broken symmetry variable u giv- 
ing displacements of the layers along z. This model also 
describes Rayleigh-Benard stripes [16], where the modu- 
lated field is the local temperature, which is not a con- 
served quantity. This case was dealt with briefly in [14]. 
The hydrodynamic, long-wavelength model for the dy- 
namics of the it-field, retaining terms permitted by sym- 
metry, including u — ¥ —u, z — > ~z [26], to leading order 
in gradients and in powers of u, reads 

dtu = Bd 2 z u + DV\u - KV\u + f u , (1) 

where /" is a Gaussian, zero-mean spatiotemporally 
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white noise with variance 2A. The term with coeffi- 
cient D [27] is forbidden by rotation-invariance of the 
free energy in an equilibrium smectic without an align- 
ing field. It is however permitted here simply because 
rotation-invariance at the level of the equation of motion, 
which is all one can demand in an active system, does not 
rule it out. Its physical content is that layer curvature 
produces a local vectorial asymmetry which must cause 
directed motion of the layers as this is a driven system. 
Symmetry does not fix the sign of D. An undulation 
instability [28] arises if D < 0. For positive D the lin- 
earity and spatial homogeneity of (1) makes it straight- 
forward to show via spatial Fourier transform that the 
variance (|u(q, t)| 2 ) = A/(Bq 2 + Dq\) in Fourier space 
for small wavevectors q = (q_i_,<7 z ). This implies [14] 
a real-space variance (u 2 ) = J" (|u(q, t)| 2 ) that is finite 
in d — 3, corresponding to long-range smectic order; in 
(1 = 2, (u 2 ) ~ logL for system size L, corresponding to 
quasi-long-range order. This establishes result (i) for the 
simplest case. 

Ignoring the K term, which is irrelevant at long dis- 
tances, rescaling rj_ = r'^D/B) 1 / 2 , z' = z, and express- 
ing u in terms of the angle 6(r') = (2ir / a)u(r) where a is 
the layer spacing, we can rewrite (1) in the form 

d t 9 = BV 2 6 + f s , (2) 

with rescaled noise statistics 

d-l 

(f(r',t)f(0,0)) = ^A^j 2 5(r' x )5(z')5(t). 

(3) 

Defining k = B^-^^D^^^a 2 /(2tt 2 A) it can be 
shown that the steady-state probability distribution for 
9 implied by (2), (3), is exp[(-/«/2) / dV(V'<9) 2 ], identi- 
cal to that for a thermal equilibrium XY model with a 
stiffness/temperature ratio k. 

This equivalence to an equilibrium XY model implies 
that topological defects (i.e., dislocations) in an active 
smectic in dimension d = 2 unbind, driving the system 
into the active nematic phase, [20] when k = 2/tt, i.e., 
when 2?r 2 A/ a 2 (BD) 1 ' ' 2 = vr/2 . This locus is plotted in 
the A-D plane in figure 1(a). 

Since it is a purely active effect, we expect D oc cq. 
We also expect the noise strength A to get an active con- 
tribution proportional to Co, and a co-independent ther- 
mal contribution proportional to ksT. Hence, varying 
Co maps out a straight line with positive intercept on the 
A-axis in the A-D plane, as illustrated in figure 1(a). 
As is clear from that figure, this experimental locus can 
only enter the active smectic region by crossing the active 
smectic to active nematic phase boundary twice. Hence 
our conclusion that re-entrance is inevitable in two di- 
mensions for these systems. 

In three dimensions, the situation is quite different, be- 
cause equilibrium smectics are stable against dislocations 




(a) d = 2 (b) d = 3 



FIG. 1. Phase diagram of active smectics, in (a) d = 2 and 
(b) d = 3. 

in d — 3, since dislocations are line defects that remain 
bound even in the absence of the D term in (1) [29] As 
a result, the active smectic- active nematic phase bound- 
ary does not go all the way down to A = at D = in 
d = 3. It does, however, develop an infinite downward 
slope at D = 0, as can be deduced by the following argu- 
ment: for small D, as the transition is approached, the 
system will act like an equilibrium (D = 0) system until 
the D term in (1) becomes comparable to the K term at 
in-plane wave vectors q± ~ £j j where £j_ is the equilib- 
rium, in-plane correlation length for smectic order. This 
condition leads to -D/£± ~ K/£±i which implies 

U ~ y/I/D. (4) 

Near the equilibrium AN transition, £j_ oc |T — 
Tan\~ u± where v\_ is the "thermodynamic" equilibrium 
correlation length exponent in the _L direction; see ref- 
erence [30] for a further discussion. Since A plays the 
role of temperature here, (4) leads to a shift oc D x l 2vi - 
in the critical A cr . Both theory [30] and experiment [31] 
find 1/2^_l < 1, so the phase boundary in figure 1(b) has 
infinite slope as D —> 0. 

The locus in the A-D plane mapped out by varying 
Co remains a straight line; now, however, re-entrance, 
though still obviously a possibility (e.g., for locus 2 in 
figure 1(b)), can be avoided, as on locus 3 in figure 1(b). 

As in d — 2, in d = 3 the active smectic to active ne- 
matic transition is in the XY model universality class for 
the model with no conservation laws. While the other 
four dynamical models we study in this paper are not 
equivalent to any equilibrium XY model, they have ex- 
actly the same scalings of their equal time u-u correla- 
tions as this simplest model. Hence, we expect similar 
phase diagrams, with the critical A cr oc D 1 ! 2 for small 
D in d = 2. However, the universality class of the phase 
transition in these other models may be different. 

We now turn to our second model, in which the number 
of active particles is conserved, but momentum is not. 
Now the concentration c of active particles with mean Co 
and Sc= c — c , joins the broken symmetry variable u as 
a hydrodynamic field. For apolar phases, the equations 
of motion must be unchanged if u — > —u and z — > — z 
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simultaneously. The c-dependent term in the u equation 
of motion that is lowest order in spatial derivatives that 
respects this symmetry is oc d z c; hence, the equation of 
motion for u is: 

d t u = Bd 2 u + DV\u - kv\u + Cd z 8c + f u , (5) 

with /" having the same statistics as those in (1). 

The equation of motion for the concentration c must, 
due to conservation of total number of particles, be ex- 
pressible in the form: dtc = — VJ C . Gradient expanding 
the current subject to the symmetry constraints gives: 

J c = -z[(A z - W c )d z Sc + W c c W 2 ± u + C zz d 2 z u] 
-V ± [A±Sc+(C ±z + W c c )d z u} + f c (6) 

where the Gaussian noise f c = (f±,f z ) has statistics 

</f(M)/?(O,0)) = (Ai^ + A c z 5* 3 )6 d (r)5(t). (7) 

In (6) we have included an active current [5] W C V • (cnn) 
with n = (z — Vu)/|z — Vti|, where W c is a phe- 
nomenological coefficient. In an equilibrium two com- 
ponent smectic, (6) and (7) would hold, but with the 
constraints W c = and C ±z /C zz = A ± /A z = A c ± /A c z . 

The two key results that emerge from (5)-(7) are that, 
for q — > and for all directions of q, the equal-time cor- 
relators (|w q | 2 ) oc 1/q 2 for all directions of q, and that 
(|<5c q | 2 ) is finite. As before, this q~ 2 scaling of u fluctu- 
ations implies translational order is quasi-long-ranged in 
d = 2, and long-ranged in d = 3. The finite concentra- 
tion fluctuations, which result from the smectic elasticity 
B, imply the absence of giant number fluctuations (our 
result (iv)). 

We next consider active smectics suspended in a fluid 
medium. The total momentum of suspended particles 
and ambient fluid is conserved; the corresponding mo- 
mentum density g is therefore slow and hydrodynamic. 
The other hydrodynamic fields u and c remain as well, 
of course. We'll assume overall incompressibility, so that 
total (particle + fluid) mass density p = p = constant 
and V • v = 0, where v = g/p is the velocity field. 

Conservation of total momentum reads d t g = — V • cr, 
with a linearized stress tensor 

er = pi - »7(Vv + Vv T ) + er (e/ > + <r a + <t n , (8) 

with p the fluid pressure, rj the viscosity tensor, and the 
elastic force density —V ■ cr^ ~ —nSF/8u, with 

F=^Jd d x [B(d z uf + K{V\u) 2 + A(5c) 2 + 2C6cd z u] . 

(9) 

Here B and K are layer compression and bend mod- 
uli respectively, A the osmotic modulus, and C a cross- 
coupling. The active stress [3, 8] <r a = — Wain, with 
negative and positive W corresponding respectively to 
extensile and contractile stresses, while <t n is noise. That 



<r a ex c follows because W is the activity per particle. As 
in any smectic A, the layer normal n is geometrically 
locked to the displacement field: <5n = n z ~ — Vj_it. 

The resulting equation of motion for v, linearized in v, 
u and 5c = c — Co, with Co the mean concentration, reads 

Po d t v = -Vp + i[Bd 2 u - KV\u + (C + W)d z 5c] 

-Wc a {iV 2 L u + d z V ±_u) -r-v + r, (10) 

where — T v = V • (r)Vv), f™ = V • cr N is a 

momentum-conserving noise, and (<7^(0, 0)cr^(r, <)) = 
2AijkiS(r)S(t). Aijki is uniaxial, symmetric in ij and kl 
and under interchange of ij with kl, and thus has five 
independent components in d — 3. In thermal equilib- 
rium Aijki oc r/ijki, but not in general nonequilibrium 
systems [32]. We take A ijU = A v (6 ik 5ji + 5 a S jk ) and 
Vijki = V^jk^u- The linearized hydrodynamic equation of 
motion for u is 

d t u = v z + Bd 2 u + DV\u - KV\u + Cd z Sc + /"; 

(11) 

where the noise /" has statistics as in (1) [33]. 

The equation of motion for cis: d t c = — V • J c with J c 
given by equation (6). 

For simplicity, we take A c ± = A L Z = A c , C ±z = CA C = 
C zz , and A± = AA C = A z , where A and C are as in 
(9), and define CA C = E, AA C = G. Activity now enters 
only through the W terms in (10). 

To study the Stokesian limit, we neglect inertia and 
acceleration, impose incompressibility V • v = in (10), 
and solve for v in terms of u and c. Inserting the result 
in (11), and defining $ = —d z u, we find that the spatial 
Fourier transforms <5c q , <f> q obey 

9 t $ q = -M^{[Bq 2 z + Wc {q 2 z - ql) + A^i]$ q 
-{C + W)q 2 5 Cci } + [5 t $ q ] P - iq z (f% + Mq/y (12) 

d t <Sc q = (Eq 2 + 2VF c c <zi)$ q - (Gq 2 - W c ql)5c^ -iq f^ 

where M q = q\jr\q^, and [<9 t <I> q ]p summarizes the "per- 
meative" B, K, C , D terms from (11), which are of higher 
order in wavenumber than those shown explicitly in (12). 

Suppose B > CE/D = C 2 j A, so that when activity 
W = the smectic state is stable. Let |W| > C > 0; 
a similar analysis holds for C < 0. At small q, where 
[dt^ci\p is negligible, it is clear from (12) that negative 
(i.e., extensile) W, can lead to an instability with q along 
z, i.e., a modulation in layer spacing. However, the layer 
compression modulus B always stabilizes this when (B — 
\W\cq) > 0. Thus, the system is stable for small enough 
\W\, establishing our result (iii). 

For contractile active stresses W > 0, we see from (12) 
and (13) that the most unstable modes have q in the 
_L direction, in which neither the layer compression elas- 
ticity nor the coupling to the concentration act. Hence, 
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the instability threshold for W — > in the limit of large 
system size, as in [8, 9]. The instability causes splay and 
self-generated flow, as in active nematics [8, 9]. For smec- 
tics, this is a spontaneous version of the Helfrich-Hurault 
[15, 28] undulation instability. 

We turn next to the effects of confinement. Consider 
an active smectic with layers normal to the z direction, 
confined between no-slip walls parallel to the xz plane, 
a distance I apart. We start with (10) but with both P 
and r ~ 77 /£ 2 = T nonzero at zero wavenumber because 
the walls are a momentum sink. Solving the modified 
(10) for the (now) fast variable v in terms of the slow u 
and c and inserting the result in (11) yields 

a t $ q = - (Bq 2 + Dq 2 x + Kq 4 x )^ + Cq 2 <5c q 

+Ar q £ q ), (14) 

while equation (13) continues to hold. Here B = B + 
(B + WcojM^, D = D- WcoMq, K = K + KMq, and 
C = C + (C + W)Mq. The wave vector q now lies in 
the xz plane and the mobility M q = (ql/q 2 )(l/T), in 
contrast to the bulk system, does not diverge at small 
q, so that terms involving M q , which arise from the 
(screened) hydrodynamic interaction, are of the same or- 
der in wavenumber as the permeative terms [<9t < f> q ]p. Ig- 
noring the concentration field, we see that the relaxation 
rate of layer displacements with wavevector in the x di- 
rection is now proportional to [D — (c$W /T)q 2 / q 2 ]q 2 for 
q x > q z and {B + [(B + c W) /T]q 2 x / 'q 2 }q 2 z for q z > q x . 
Thus, there is a range of parameters for which the active 
smectic is stable (result (vi)). For other parameter ranges 
instabilities occur; e.g., if D — in (11), an undulation 
instability occurs for q _L z, despite confinement. 

The instability that arises in the extensile (W < 0) case 
when |W| > B is interesting. Equations (14) and (13) 
have the same form as the linear part of the Fitzhugh- 
Nagumo [34, 35] equation, which exhibits sustained os- 
cillations under rather general conditions. We speculate 
that such oscillations could also occur here; i.e., a breath- 
ing smectic. We will explore this in future work. 

We now turn to fluctuations in the bulk Stokesian 
limit. For q z — > 0, 8c drops out of (12), and it is 
then easy to check for extensile activity that (|u q | 2 ) = 
2A t ,/|IU|cog 2 ^. Somewhat more tedious algebra, which 
we'll present elsewhere, shows that (|u q | 2 ) ~ 1/q 2 for 
all directions of q. We can also see from (12) that, in 
the Stokesian regime, time correlations of u decay at a 
nonzero rate in all directions for q —> 0. These facts 
imply that, in bulk active smectic suspensions as well, 
the coupling to u via the E term in (13) won't lead to 
diverging concentration fluctuations in general. (|<5c q | 2 ) 
does diverge, however, upon approaching the extensile 
instability, with a correlation length ~ (B — |VF|co) -1 / 2 . 

We conclude with second sound. Leaving the steady 
Stokesian regime, taking acceleration (9*v) into account, 
and working at long wavelengths where viscosity is neg- 



ligible, makes the coupled dynamics of v and u, for 
wavevectors in the plane of the layers poftv = — Vp — 
JUcozV^u; d t u = v z ; V • v = 0, can readily be seen to 
give sound waves with a speed \J— Wcq/pq. 

In conclusion, we have constructed the dynamical 
equations for active smectics, both in bulk suspensions 
and in confined systems in contact with a momentum 
sink. Our theory is generic, applicable to any driven 
system with spontaneous stripe order. We show, extend- 
ing [14], that noisy active smectic order is long- ranged 
in dimension d = 3 and quasi-long-ranged in d = 2 
for all dynamical regimes, and that active smectic sus- 
pensions have a nonzero second sound speed parallel to 
the layers. For d = 2 we predict a Kosterlitz-Thouless 
transition from active nematic to active smectic, with 
a re-entrant nematic at low concentration. We show 
that smectic elasticity suppresses the giant number fluc- 
tuations and extensile instabilities that occur in active 
nematics, but that bulk contractile systems exhibit an 
active undulation instability. Active extensile stresses, 
if strong enough, give rise to a "breathing" instability 
which is likely to be oscillatory. Our results should ap- 
ply to a wide range of active systems, including horizon- 
tal layers of granular matter agitated vertically or fluids 
heated from below. 

We look forward to detailed experimental tests of our 
predictions. 
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